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Abstract 

In this report we present the detail theory of the focusing in bent crystals using the 
statistical mechanics. The quantum mechanical effects for focusing and abberations 
are not considered. The beam structure in phase space are consedered and envelope 
near the focusing spot and intensity profile are derived. 
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1 Introduction 



When a high energy particle with a momentum P are captured in a channel- 
ing motion of bent crystal, the equivalent magnetic field acting on the particle 
can be estimated as ~ P/(0.3R), where R is a curvature radius of the bent 
crystal [1,2]. For the momentums P > 300GeV/c and radius R ~ lm the 
equivalent magnetic field can reach the values > 1000 tesla. This phenomena 
has been extensively investigated for a beam deflection, splitting, extraction, 
and spin precession measurements [2, 3, 4]. The first experiment of the bent 
channeling[5] confirmed that the bent crystal could also be used for construct- 
ing a focusing element with extremely short focal length and small focal spot. 
Several crystal devices for focusing were suggested[6, 7,8,9], but the successful 
experiments[10,ll] were done with the crystal having a special shape drawn 
in Fig. 1(a). This is a regular bent crystal with the output face carved in such 
a way that the tangent lines to the planar channels on that face are converged 
to a focal point. With these conditions, the center points of the channels on 
the output face of the crystal constitute a cylindrical surface with diameter 
D, greater than curvature radii of bending, D > R ma x > R > Rmin- The 
estimations [12] show that the size of smallest focal spot is proportional to the 
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Fig. 1. Geometry of focusing crystal for high energy particles (angles, thick- 
nesses and curvatures are exaggerated), (a) The radius-vectors Rq, Ri, ...Rn orig- 
inated from center O designate the curvature of each plane. (b)The n-th channel 
(n = 1, 2, 3...N) is located between n-th and (n+l)-th planes and has a 'local' Carte- 
sian coordinate system (x', z f ) with offset z n ) and angle direction # n = <p n — <p Q 
relatively to 'base' Cartesian coordinate system (x, z) of the channel 0. 

square root of the crystal thickness d c . There is also an evaluation of the in- 
tensity at the focal spot, and the crystal geometry is proposed where the focal 
parameters can reach the extreme values. 

In this article we present the detail theory of the focusing in bent crystals 
using the statistical mechanics. The quantum mechanical effects for focusing 
and aberration are not considered here. 



2 The model potential of focusing by bent crystal 

We consider the classical model of channeling [13] with a continuous tem- 
perature dependent planar potential, when ions are captured into channel- 
ing trajectories and directed along the bent atomic channels. Similar consid- 
erations in bent crystals without boundary conditions were done by many 
authors [14,15,16,17,4] and based on the assumption that curvatures of crys- 
talline structure are small and bent plane locally looks like flat surface. In local 
coordinate system the model of continuous straight planar potential [13] can 
be apply and describe the most important features of bent channeling [4] . If 
each plane has a constant curvature, it is convenient to use a polar coordinate 
system (p, if) with a center located at the point O (see Fig. 1(a)) and direction 
(p, ip = 0) tangent to the cylindrical surface. In such coordinate system the 
crystal potential for the bent channeling in Moliere's approximation [18,19,20] 
with a different length of the channels can be described as following 



N 



U b -c(p, V 3 ) = U n(P ~ Rn, <p), 




n=0 
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where R n = R m in + n * d p is the radius of curvature of n-th bent plane (Rq = 
R m i n ,RN = Rmax), d p is interplanar distance, U n (p — R n ,<p) is temperature- 
dependent continuum potential of this plane derived from the Moliere poten- 
tial [21,22] for straight plane: 



(exp(A5)er/c(A^ + ^) 

+ exp(-/?4)er/c(/3 i ^ - ^)), ip min < ip < ip n , 

0, <p> < (fmin, 



(2) 



Here a and w are Thomas-Fermi screen radius and root-mean-square tem- 
perature displacement of crystal atoms (in Si at 293K u = 0.075A), a; = 
{.1, .55, .35}, Pi = {6.0; 1.2; .3}, U = 7iZe 2 an p (n p is a density of atoms in 
the plane), erfc(x) is the complementary error function. The angles <p m i n , 
denote the enter and exit surfaces of the nth bent plane. Thus, the end cylin- 
drical surface of the crystal lies in the range ipo < ip < (p ma x- Others notations 
can be found in Fig. 1. Although Moliere's approximation (2) for isolated 
plane has a simple analytical form, it can not directly be implemented for 
calculation of the channeling trajectory, because the channel potential is su- 
perposed from several neighboring planes. However, if we consider an infinite 
crystal at fixed temperature, this superposition can be fitted by a polyno- 
mial of suitable order. The estimations can be easily received by an harmonic 
approximation [4], which coincides with Moliere's approximation on the crys- 
talline planes and on centers of the channels. By involving 10 crystalline planes 
located far away from boundaries (K » 1), the parabolic potential can be 
fitted to the Moliere's approximation with the following parameters: 

K+5 

U m ax= U n ((K -n)*d p ,ip), 

n=K-4 

K+5 

U min = U n{{K ~ Tl) * d p - dp/2, ip) . (3) 

n=K-A 



One can extend this model to the crystalline potential near the boundaries 
ignoring the fact that potential of the first several subsurface channels becomes 
slightly asymmetric (see Fig. 2). Outside the crystal, the potential can also 
be fitted by half part of parabola, for example, with a condition that on the 
boundary plane the potential is continuous and the parabola intersects the 
Moliere's potential at the value U m i n /2. Because other planes give a small 
contribution to the surface potential, the distance of this intersection from 
boundary plane is very close to d p /2. The resulting system of two equations, 
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Fig. 2. The crystal potential for the range <p m i n < <p < *£o in harmonic approxima- 
tion with parameters U max = 35.9eV, U m in = 13.3el/, and Moliere potential for 
Si (110) at T=300K. Notations: q = (p - R m i n )/d p , U nrm = U max - U min . 

U b * ( Pb /d p + 1/2) 2 = U min /2, 

Ub * (Pb/dp) 2 = U max , 

leads to the boundary semiparabolic potential with parameters: 
(2U„ 



(4) 



Ub 



-U ■ ) 2 

' max v mm I 



U ma x ~T" U m i n /2, -\- \J 2 U max U n 



Pb = -d p 



Umax ~T" yUmaxUmin/2 

2 U max Umtp 



(5) 



Finally, the parabolic model of Moliere's potential (2) of the bent crystal with 
boundaries can be written (see Fig. 2) 



r ^ . U b [ P R l Pb ) , p b + R <p<Ro, Vmin < ^ < <P0, ^ 

0, otherwise, 



- for boundary and 



Un(p, <p) 



TT ( 2 (p-Rn- d p/ 2 ) \ 2 I J J . p ^ _ ^ p 

nrm ^ j J t y-> mini ±L n ^ y ^ -"-n+lj 

Pmin < V < <Pn, 

Ub ^ P -R n ^- Pb y ) pb + Rn+i <p< Rn+u 

ip n < (f < <f n +l, 

0, otherwise, 



(7) 
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- for channels < n < N — 1. The boundary potential (6) does not play any 
role for bulk channeling, but may play an essential role for an abberation of the 
focusing, if the length of this potential along the n-th microbeam, R n (ip n+1 — 
(p n ), greater then the length of oscillation of particles inside the channel, A c ~ 
d p /typ is the Lindhard critical angle for planar channeling). This condition 
may happen for extreme focusing with shortest focal length [12] (see also 
below) and a similar potential term at the end of each channel must be kept 
for corrections of phase curves (see second line of (7)). If this condition does 
not hold and distortion caused by boundaries is negligible, and the channel 
potential (6,7) becomes: 



N-l 



Ub-c(p,<f) = J2 U n(P,tp), (8) 



n=0 



where 



U n (p,ip) = < 



TT ( 2 (p- R n~ d p/ 2 ) \ 2 _J_ T J D ^ n - D 

u nrm y ^ J t ^rnini lh n ^ \ J ^ ll n+l: 

n ■ (9) 

0, otherwise. 



3 Dynamic of particles inside the bent crystal 



Consider the motion of a relativistic spinless particles in potential (1) or (8). 
The Hamiltonian of the particle is defined by 

H(Q, P) = cy/PiPi + mW + U b . c (p, </?), (10) 



where Q = {p, y?} and P = {p p ,Pip} = {mp / -f,mp 2 Lp'-f} are canonical coor- 
dinates and conjugate momenta of particle, 7 = (1 — (p 2 + p 2 (p 2 )/c 2 )~ 1 >' 2 , 
PiP 1 = p^+p^/p 2 - The Hamiltonian H(Q, P) is a piecewise function of ip and p 
and, most important, does not depend on ip inside and outside of each individ- 
ual channel. Hence, the angular momentum of particle p^ = mp 2 tp^ (see Fig. 1) 
is conserved inside and outside the crystal and may change only when particle 
crosses the boundaries of the crystal. The phase curves of the particles are de- 
scribed by system of Hamiltonian's equations ( P = —dH(Q,P)/dQ, Q = 
dH(Q,P)/dPy. 
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Fig. 3. The effective crystal potential for bounded crystal Si (110) at T=300K in 
harmonic approximation, Eq. (6,7), for radial sections: (p m i n < cp < <po - solid line; 
(fa < ip < cpi and cp\ < cp < cp<± - dash lines. In the bent crystal, in addition to 
the normal channels, particle may also have a surface channel in which it moves 
sequentially reflecting from a surface plane. Note that this surface channel is much 
wider than the crystal channel and might play an important role for focusing in 
extreme conditions [12]. 

• c Pp • c Vip (11) 

P ~ VP t P l + m?c? ' V ~ p 2 ^/P~P l + m 2 c 2 ' 1 ' 

Using the fact, that last two equations of (11) can only be resolved if 

c 1 , „ 

(12) 



\/PiP i + m 2 c 2 Twy 
the system (11) can be written in the form: 

P 2 

Pp = — - dU b - c (p, <p)/dp, = const, 
m^p 6 

P=^ <P=-^- (13) 

The Hamiltonian (10) of the system (13) has a cylindrical symmetry, therefore 
the angular momentum of the particle p v is conserved. The tangent momentum 
of particle p^j p along the curved channel is not exactly conserved, but taking 
into account that the channeling particle can not move further than d p /2 from 
the center line of the channel and d p /R n -C 1, the tangent momentum of 
particles p t = p^/p = mptp'j and tangent velocity v t = pip = pt/{m~y) are 
conserved with high accuracy. We also make the following 

approximations: 

p p < (pt,mc), 



6 



m7p 3 ?niR m 



U b _ c (p, <p) < c^/p^ + m 2 ^, (14) 



which mean that the incident particles are directed at the small angle to the 
crystalline planes, the centrifugal force p^/ijn^p^) = p t v t /p for all crystal 
channels is practically the same (the crystal thickness d c = Nd p <C R m in), 
and gamma factor inside the crystal is changed by a small amount: 



7 - 7o = 



\jv 2 po+VlJpl + m2c2 



mc 



(15) 



From (11), (14) and (15) the well known equations [1] of motion of a particle 
in bent crystal can easily be received : 



p 2 

m>y p = |g dU b ^ c (p, (p) /dp, = const, 

m 1 'o-K-min 



m lo R m in m lo 



Note that parameter p in Eq. (11,13,16) means only the piecewise character 
of potential (1, 8), which has discontinuity at the front and end surfaces of 
the each channel. The normal derivative of the potential dUb- c (p,(p) /dp on 
these surfaces does not exist. Of course, the real atomic potential is smooth 
function on the surface. Its derivative rendering the normal force causes an 
impulse in normal linear momentum when the particles cross these surfaces. 
The linear momentum tangent to the front and end surfaces of the channels 
does not change as particle passes the boundary of the crystal, but angular 
momentum may slightly change due to force acting perpendicular to these 
surfaces. Thus, the boundary conditions for Eq.( 11,13-16) on front surface of 
the crystal, <p = ip m in, can be presented in the form 

Po = Pi, Ppo=Ppi, Ptpo = P<pi + A P , (17) 



where {p ,Ppo,Pipo), {PiiPpiiP^i) denote the radial coordinate, linear tangent 
and angular momenta on the particle for p — > y?mm and tp — > ip min respectively. 
Using the law of conservation of the energy 

c^Jp 2 o + Vlo/Po 2 + m 2 c 2 = c^p 2 pi + pli/pi 2 + m 2 c 2 + U h - C (p u p min ), (18) 
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we can receive a change in angular momentum p^ in (17), 
^ PyP 2 po + Plo/p 2 o + m 2 c 2 

A p = i U b - C (pi, (Pmin), (19) 

on the boundary of the crystal. The same boundary conditions (17) take place 
when the particle passes the end surface of the crystal. The assumption (14) 
allows also to split the energy of relativistic particles into kinetic energy of 
particle along the bent channel E\\ and transverse energy E± [23,24,14,25,17]: 



E = En + E A 



The transverse energy E± is constructed of three terms: the kinetic transverse 
energy, the static potential energy and centrifugal potential energy. The last 
two terms may be combined in one - the effective transverse potential in a 
bent crystal (see Fig. 3)): 

UellM = UM -^z^l. (21) 



Since the energy along the bent channel Eu is conserved, the transverse energy 
Ej_ is also conserved. Equation (16) with the potential (8,9) and boundary 
conditions (17,19) can be easily resolved, 



p= {po- Pc,n) COs(^t) H P — Sin(atf ) + p C) n, 

P P = P P o cos(wt) - (p - p c ^ n )ujm^ sm(ut), 

(p = Qt + ip m i n , 

p v = const, (22) 



where the frequency of angular rotation, frequency of transverse oscillation 
and the radius of the central line of motion for n-channel are 



Q = 



Py 



ITL'-foR 



■2 ' 
min 



UJ = 



SUp 



p 2 

Ftp 



d p 

Pc,n = Rn + y + 7^2^ 



(23) 
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Thus, the motion of all particles in phase space can be described by Eqs. (22) 
with parameters (23). The maximum amplitude of oscillation and transverse 
momentum inside the bent crystal are (see Eq.(22) and Fig. 3) 



| Prnax Pc,n \ 



dp 



V„ = (f - ^ m2 * lRl )*m l0 . (24) 



^mtn 



To facilitate further calculations we introduce the normalized dimensionless 
local coordinates x, v x for each channel: 



dp/ 2 oom^odpl 2 



which satisfy the equation of conservation transverse energy (20) 

x 2 + v 2 x = r 2 ch . (26) 
The parameter r c h, 



^ 2 m 2 ^R min d p e c ' 
V% e = 4E/q 

m loR m in dp 



e= — d3 > e c = ( 27 ) 



can be called the dimensionless radius of the microbeam in the phase space, 
and e, e c are an effective electric field and critical electric field produced by 
bent crystal [4,3]. 



4 Statistical model of microbeam focusing 



The classical approach to the kinetics of beams assumes the Boltzmann equa- 
tion for particle distribution function / = /("f > , J?) in a steady state situation: 

lf % + e ^% = St (f)> ( 28 ) 
or op 
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where el? is electric field acting upon the ion with charge e inside the crystal, 
St(f) is a collision integral, which plays an important role in establishing 
the equilibrium states and dechanneling inside the crystal. Outside the crystal 
St(f) and e~E vanish and the Boltzmann equation reduces to particular simple 
form 

= °" ( 29 ) 



Consider 2-D geometry for plane channeling (Fig. (la)). The particle beam has 
an initial distribution fi n (x, z in (x),v x , v z ) on the entrance of the crystal surface 
z = z in (x) and some distribution f ou t(x, z out (x), v x , v z ) formed by channeling 
process on the end face of the crystal z = z out (x). We use Cartesian coordinate 
system instead of the cylindrical coordinates used in Sec. 2, 3 and assume that 
all coordinates and their conjugate momentums or velocities are dimensionless. 
Behind the crystal the Eq. (29) becomes 

df(x,z,v x ,v z ) df(x,z,v x ,v z ) 
v* g~ x + v, Wz = (30) 



with general solution 

f(x, z, v x , v z ) = g(x- z* v x /v z , v x , v z ), (31) 



where g is arbitrary function of three variables. The boundary condition on the 
end face of the crystal z = z out (x) gives unique solution for beam distribution 
(v z = const): 

g(x - z out (x) *v x /v z ,v x ,v z ) = fout(x,z out (x),v x ,v z ). (32) 



As a simple model, consider single channel, say channel with output angle 
ipo (see Fig. la). Output particle distribution on the end face of the channel 

(zout = 0) is 



fout,o(x,0,v x ,v z ) = h(v z 



, v 2 x + x 2 < r 2 ch , 

ch 

0, v 2 x + x 2 >r 2 ch , 



(33) 



where h(v z ) is some distribution in space of longitudinal component of the 
velocity, r c h is the radius of particle distribution in the phase space which is 
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Fig. 4. Evolution of distribution function for microbeam from single channel. The 
trajectories of particle in 2D are parallel to the x axis. The phase profile along 
trajectories of particles moves with a constant velocity (v x = const). 

less than d p /2 radius of crystal channel. The distribution (33) is normalized 

2 

dp OO OO 

J fout,o(x, 0,v x ,v z )dxdv x dv z = 1, (34) 

2_ — OO — oo 

dp 




viz., one particle per one channel. By using boundary conditions (32) and (33) 
it is easy to show that the distribution function of this channel in phase space 
behind the crystal will be 



fout,o(x,z,v x ,v z ) = h(v z ) < 



^,v 2 x + (x- z v ff < r 2 ch , 
0, vl + (x - z^f > r* ch , 



(35) 



and the evolution of this distribution function along the axis z is shown in 
Fig. 4. In phase space, the cross section of distribution is deformed, but the 
area of microbeam is conserved as well as the phase density along trajectories 
in the agreement with Liouville's theorem. 

The same treatment can easily be extended to each channel if we take proper 
Cartesian coordinate system (x, z) with origin at the center of this channel on 
the end face of the crystal and with OX axis directed to the focus (see Fig. lb). 
The coordinate system for n channel is connected to Cartesian coordinates 
(x, z) of the 1st channel by following relations 



x -x n = x n cos($ n ) - z n sin($ n ), 

Z - Z n = X n sin($ n ) + ZnCOsidn), 



(36) 



where the origins of (x, z) coordinate systems are 
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Fig. 5. The phase space of microbeams at focal spot when F m i n « F max . Only ten 
channels are taken for illustrative purposes. In fact, about 10 7 crystal (110) channels 
would typically be involve in forming a focal spot for the Si crystal with thickness 
d c = 2mm. 

X n {Rmin Tldp -\- dpl'2^COs( y '& n ) Rmin) i 

z n = (Rmin + nd p + d p /2)sin{d n ), (37) 
and the angles i? n = (p n — (p denned by (see Fig. la) 



{Rmin ~\~ Tidp) y-D 2 Rmin Rminy D 2 {Rmin 4" ndp) 2 

sin{V n ) = — 1 . (3£ 



Among others important parameters of the focusing (Fig.l), there are simple 
formulas for the maximum and minimum focal distances 



D 2 - R 2 

- 1 "mm > 



Rmin \J D {R mm -\- d c ) , 

d c = Nd p . (39) 

The distribution function for n-th channel in phase space behind the crystal 
is 



1 n|2 I ( ry. 4 V± n \i ^ r 2 

fout,n{ X n,Z nj V ±n ,V £ J = h(v £ J < ^ ^ (40) 



0, vj +{x n -z n v ^) 2 >rl 



ch-> 

2n 



The equations (36) represent a reversible transformation, i.e. that they also 
define the {x n , z n ) as functions of the (x, z), or, in other words, that they are 
soluble with respect to the (x n , z n ) 
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x n = {x- x n )cos{& n ) + (z - z n )sin($ n ), 
z n = -(x - x„)sffl(i) + (z- z n )cos{-& n ). 



(41) 



For the distribution function (35) we have also to calculate v Xn and V£ n which 
can be done by differentiating (41) with respect to t: 



= v x cos($ n ) + v z sin('d n ), 
-v x sin(§ n ) + v z cos(§ n ). 



(42) 



Note that the distribution in longitudinal space of the velocities is invariant 
under the changes of (42), i.e. h(v z ) = h(v Zn ). Substituting in the equation 
(40) the values of the (x n , z n , v Xn1 v Zn ) given by formulas (41), (42), we get the 
equation for distribution functions for all channels in coordinates (x, z,v x ,v z ) 



fout,n{x,z,v x ,v z ) = h(v z ) < 



i 

■rrr' 2 



(v x cos(ti n ) + v z sin($ n )) 2 + 

{{x - x n )cos($ n ) + (z- z n )sin($ n )- 

((x - x n )sin(§ n ) -(z- z n )cos(§ n )) 

)+v z sin(-d n ) 
v x sin($n)-v z cos (i?„ ) 



< r 



chi 



(43) 



(v x cos(ti n ) +v z sin{d n )) 2 + 

[{x - x n )cos($ n ) + (z- z n )sin($ n )- 

((x - x n )sin(d n ) -(z- z n )cos(§ n )) 

v x cos(tf n )+v z sin(tf n ) \ 2 y r 2 
v x sin($„)—v z cos(i)„)) — chi 



The final distribution function for the beam behind the crystal is 

N-l 

fout(x, Z, V x , V z ) = fout,n( X i Z > V x, V x ), (44) 

n=0 



where N = d c /d p - number of plane channels of crystal. If we substitute 
z — F max in (44), we can get the distribution function in a focal plane of the 
phase space. The main features of this distribution can be seen in Fig. 5. All 
ellipse have the same area as mentioned above. The upper ellipse is shaped by 
microbeam from channel 0. It has greater x-size because the larger distance 
from the output channel to the focal point. The bottom ellipse is related to 
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the last microbeam. It has a smallest deformation. The centers of the ellipses 
in focal plane are located on the V^-axis at the points 



v z sin($ n ) 
cos(d n ) 



(45) 



For an experiment measurements, it is important to find the intensity profile 
at the focal spot. It can be received by the integrating Eq.(44) over v x : 

fouti^i F max ) / / y ' fout,n(%i F max ,v x ,v z )dv x dv z . (46) 

J J n=0 



In practical calculations, the summation over a pile of channels (see Fig.5)can 
be change by integrating over n.To simplify the calculations, we assume that 
all i? n << 1. For normalization (34) this gives the pick value of distribution 
function in the form 



fouti^i Fmax) — 2 r c y t V 



v z 2 + D 2 cos(T) 2 



D d p sin (r) 

2 + D 2 cos(r) 2 



D d p sin (r) yjv z 2 + D 2 cos (r) 2 - 2 D d p sin (r) N, 



r c h v z N , , „, 

+4 , (47) 
yjv z 2 + D 2 cos (r) 2 — 2D d p sin (r) N 



where the following parameters were introduced 



sin{Y) = — 



JD 2 - R 2 min 

cos (T) = 1 (48) 



Note that in order to calculate the envelope of the beam in phase space, we 
may take a derivative of microbeam surface with respect to n considering 
variable n as continuous and eliminate n from two equations 



$ n (x,Z,V x ,V z ) = 0, 

d$ n (x,z,v x ,v z ) 



dn 



0, 



(49) 
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where $ n (x, z, v x , v z ) = is the equation of the microbeam surface in phase 
space. We can substitute z for F max (see Eq.(39)) in equations (49), so getting 



{%) F maX) V x , Vz) 



an °" (50 » 



the equations of the cross section envelope line at the focal spot. If the distri- 
bution (40) or (43) is used as a model, we have the implicit equation of n-th 
microbeam in focal plane (x, v x ) which actually used to depict Fig. 5: 



(Bv s — Cxf 2 r ch 2 

®n{X, F max , V x , V z ) = — — —2 +V S - - 2 = 0. (51) 

(v s sm(i} n ) - C) l-sin(i} n ) 



The coefficients B, C here are given functions of the $ n : 



C 



yjl- sin{& n ) 2 
B = Dcos{T +7? n ), 

v s = v x + C sm(i?„) (52) 



The calculation of the envelope by Eq.(50) with surface (51) are straightfor- 
ward but cumbersome. In the case v s sin($ n ) << C, which is always true 
for fast particles, the envelope has a form of two deformed hyperbolas with 
asymmetry. The parametric presentation for these curves (xi,v Xi i),(x2,v Xi 2) 
are (the parameter for curves is i? n ) 



Xl Ccos(V n ) 
2 Ccos(0 n ) 



where B = B($ n ),C = C(d n ) are defined in Eqs.(52). These formulas are 
accurate for foci located far away from the crystal as well as for crystalline 
geometry provided the maximum magnification and minimum focusing size 
proposed in [12]. The envelope calculated from them is shown in Fig. 5. 
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5 Conclusion 



In conclusion, we describe some possible limitations of focusing which can 
substantially reduce the focusing effect. The smoothness of the cylindrical 
surface is of great importance, the amorphous layer on that surface must be 
as thin as possible, and all adjacent channels must constitute the staircase 
with step widths equals to the interplaner space. Of course, it is not always 
possible and deviation caused by surface roughness or mosaic blocks with 
widths about the length of channeling oscillation will produce a broadening of 
the focal spot and greatly decrease the peak intensity. Other well known source 
of aberration to be considered for focusing is negative Gaussian curvature of 
bent crystals (anticlastic effect). The asymmetry of centripetal dechanneling 
in bent channels [3,4] is a strong source of asymmetry in focal spot. There is 
also a possibility of asymmetry in coherent transitional scattering when the 
particles come out with small angles and move some short period of time in 
semichannel. For long crystal, dechanneling plays an important role, and the 
maximum possible focal distance [12] when the focal spot lies almost on the 
tip of the crystal is limited by this process. 



References 

[1] E. N. Tsyganov, Some aspects of the mechanism of a charge particle penetration 
through a monocrystal, Fermilab,TM-682 (1976) 5. 

[2] R. A. Carrigan, Developments in relativistic channeling, in: S. Chattopadhyay, 
J. McCullough, P. Dahl (Eds.), Advanced Accelerator Concepts. Plenum, New 
York, 1997, pp. 146-163. 

[3] V. M. Biryukov, Y. A. Chesnokov, V. I. Kotov, Crystal Channeling and Its 
Application at High-Energy Accelerators, Springer, Verlag Berlin Heidelberg, 
1997. 

[4] A. M. Taratin, Particle channeling in a bent crystal, Phys. Particles Nuclei 
29 (5) (1998) 437-462. 

[5] Y. N. Adishchev, et al., JETP lett. 30 (1979) 402. 

[6] R. A. Carrigan, On the possible applications of the steering of charged particles 
by bent single crystal, Fermilab FN-80/45 (1980) 46. 

[7] R. A. Carrigan, The application of channeling in bent crystals to charged 
particle beams, in: R. A. Carrigan, J. A. Ellison (Eds.), Relativistic Channeling. 
Plenum, New York, 1987, pp. 339-368. 

[8] C. R. Sun, Aplication of channeling to particle physics, in: R. A. Carrigan, J. A. 
Ellison (Eds.), Relativistic Channeling. Plenum, New York, 1987, pp. 379-397. 



16 



[9] A. Schafer, W. Greiner, Possible use of mixed crystals to focus ion beam, J. 
Phys. G: Nulc. Part. Phys. 17 (1991) L217-L221. 

[10] A. S. Denisov, O. L. Fedin, M. A. Gordeeva, et al., First results from study of 
a 70 gev proton beam being focused by bent crystal, Nucl. Instr. and Meth. B 
69 (1992) 382-384. 

[11] V. I. Baranov, et al., Results of studying 70gev proton beam focusing by 
bent crystal, in: XVth International Conference on High Energy Accelerators, 
Hamburg, Germany, July 20-24, 1992, pp. 128-130. 

[12] G. V. Kovalev, Focusing limits of high energy particles in bent crystal, Nucl. 
Inst, and Meth. B 207 (2003) 482-486. 

[13] J. Lindhard, Influence of crystal lattice on motion of energetic charged particles, 
Danske Vid. Selsk. Mat.-Fys. Medd. 34 (1965) no.14. 

[14] A. M. Taratin, Y. M. Filimonov, E. G.Vyatkin, S. A. Vorobiev, Phys. Status 
Solidi B 100 (1980) 273. 

[15] H. Kudo, Nucl. Instr. and Meth. 189 (1981) 609. 

[16] J. A. Ellison, Bending of gev particle beams by channeling in bent crystal planes, 
Nuclear Phys. B206 (1) (1982) 205-220. 

[17] J. A. Ellison, The theory of particle motion in straight and distorted crystals, 
in: R. A. Carrigan, J. A. Ellison (Eds.), Relativistic Channeling. Plenum, New 
York, 1987, pp. 79-88. 

[18] G. Moliere, Theorie der steuung schneller geladener teilhen 1, Z. Naturforsch 
2A (1947) 133-145. 

[19] B. R. Appleton, C. Erginsoy, W. M. Gibson, Phys. Rev. 161 (1967) 330. 

[20] D. S. Gemmell, Channeling and related effects in the motion of charched 
particles through crystals, Rev. Mod. Phys. 46 (1) (1974) 129-227. 

[21] J. C. Poizat, Remillieux, Phys. Rev. Lett 27 (1971) 6. 

[22] J. H. Barrett, Phys. Rev. B3 (1971) 1527. 

[23] A. P. Pathak, Phys. Rev. B13 (1976) 4648. 

[24] V. V. Kaplin, S. A. Vorob'ev, Sov.Tech.Phys.Lett. 4 (1978) 78. 

[25] A. M. Taratin, S. A. Vorobiev, Phys. Status Solidi B 107 (1981) 521. 



17 



